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1^ \ Abstract. In this talk, we review some recent advances in the theory of 

dynamic tides in close binaries. We particularly focus on the effects of resonances 
QQ ' of dynamic tides with free oscillation modes and on the role of dynamic tides in 

the comparison of theoretically predicted and observationally inferred apsidal- 
motion rates. 
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In close binaries, the tidal force exerted by each star on its companion perturbs 
Qh' the shape of the stars from spherical symmetry. Depending on whether or not 

Q . the binary orbit is circular and the component stars are rotating synchronously 

^ I with the orbital motion, the perturbations will be static (equilibrium tides) or 

^ . time dependent (dynamic tides). In the latter case, the tidal force induces forced 

nonradial oscillations in the component stars, which may be resonant with the 
components' free modes of oscillation (Cowling 1941). 

The tidal deformation of a binary component furthermore perturbs the 
spherical symmetry of the star's external gravitational field. For eccentric bina- 
ries, this perturbation gives rise to a secular shift of the line of the apses, the 
period of which depends on the orbital parameters and on the internal structure 
of the component stars. Comparison of apsidal-motion rates derived from theo- 
retical stellar models with observationally inferred apsidal-motion rates therefore 
provides an excellent tool to test current theories of stellar structure and evo- 
lution. Until recently, the comparison was usually carried out in the static-tide 
approximation to avoid the complications of a time-dependent tidal potential. 

Claret & Willems (2002) were the first to include the effects of dynamic tides 
in a comparison between theoretically predicted and observationally derived 
apsidal-motion rates for a large and homogeneous sample of eccentric eclips- 
ing binaries. The authors found that when the sample is restricted to systems 
with accurately known stellar and orbital dimensions, the comparison between 
theory and observations showed an excellent degree of agreement. However, de- 
spite this general agreement, some systems (e.g. DIHer, AS Cam, V451 Cyg) 
still show apparently anomalous apsidal-motion rates which are not fully un- 
derstood yet. Continuing efforts in improving both theory and observations are 
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therefore essential to further increase the merit of the apsidal-motion test to 
stellar structure and evolution. 

In the following sections we outline some basic aspects of the theory of 
dynamic tides in close binaries and highlight some recent advances in the com- 
parison between theory and observations. 

2. Basic Assumptions 

We consider a close binary system of stars orbiting around each other under 
the influence of their mutual gravitational force. The first star, with mass Mi 
and radius Ri, is assumed to rotate uniformly around an axis perpendicular to 
the orbital plane with a rotational angular velocity Q which is small enough to 
render the effects of the Coriolis force and the centrifugal force negligible. The 
second star, with mass M2, is treated as a point mass. 

In what follows, wc denote by a the semi-major axis, by Porb the orbital 
period, by n = 2 7r/Porb the mean motion, and by e the orbital eccentricity. The 
tidal force experienced by the star is derived from the tide-generating potential 
£rpW{'r,t) which we expand in Fourier series as 



(Polfliet & Smeyers 1990; Smcycrs, Willems, &: van Hoolst 1998). Here, et = 
{Ri / of {M2 / Ml) is a small dimensionless parameter corresponding to the ratio 
of the tidal force to the gravity at the star's equator, G is the Newtonian constant 
of gravitation, r is a time of periastron passage, the factors Ci^rn,k arc Fourier 
coefficients, the Yp(9, (p) are unnormalized spherical harmonics with respect to 
a system of spherical coordinates that is co-rotating with the star, and the ar = 
kn + mQ are forcing angular frequencies with respect to the co-rotating frame of 
reference. The Fourier coefficients Ci,,n.k depend on the orbital eccentricity and 
are proportional to {Ri/aY^^. They obey the property of symmetry Ci^-m,-k = 
ce,m,k and are equal to zero for odd values of ^-|- \m\. The coefficients c^^mfi are 
furthermore equal to zero for \m\ > I — 1. 

Due to the dependency of the Fourier coefficients Q^m.fc on the ratio Ri/a, 
the terms associated with the second-degree spherical harmonics are dominant. 
For the remainder of the paper, we therefore restrict ourselves to the tides gener- 
ated by the 1 = 1 terms. The only non-zero Fourier coefficients in the expansion 
of the tide-generating potential are then the coefficients C2,_2,fc) C2,o,fe, and C2,2,fe- 
These coefficients are independent of the ratio Rija and are solely determined 
by the orbital eccentricity e. 

For a given orbital eccentricity, the coefficients C2,m,fe(e) generally decrease 
with increasing values of k. The decrease is slower for higher orbital eccentricities 
so that the number of non-trivially contributing terms in the expansion of the 
tide-generating potential increases with increasing values of e. For a given orbital 
eccentricity, the coefficients C2,-2,fc(e) arc furthermore largest at values of k close 
to those for which kn k. Op, where rip is the orbital angular velocity at the 




X Yl'ie, (p) exp [i{aTt-kn r)] 
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periastron of the stars' relative orbit. For a more precise definition and in-depth 
discussion of the Fourier coefficients C2,m.,fe(e), we refer to Smeyers et al. (1998) 
and Willems (2003). 

3. Resonant dynamic tides 

From Eq. it follows that for each spherical harmonic Y2^{9,(j)), the tide- 
generating potential generates an infinite number of partial dynamic tides with 
forcing angular frequencies cj^ 7^ 0. When one of these forcing frequencies is close 
to the eigenfrequency of a free oscillation mode, the tidal action exerted by the 
companion is enhanced and the oscillation mode involved is resonantly excited. 
The possibility of resonances between partial dynamic tides and free oscillation 
modes was first suggested by Cowling (1941) and is particularly relevant for 
the excitation of free oscillation modes since their eigenfrequencies are most 
likely to be in the range of the forcing frequencies induced in components of 
close binaries (Willems 2003). 

Neglecting any nonadiabatic effects, the displacement field resulting from 
a second-degree resonant dynamic tide associated with the azimuthal number 
m and the A:-th harmonic in the Fourier decomposition of the tide-generating 
potential, at the lowest order of approximation in the small parameter e^, is 
given by 



(Smeyers et al. 1998). Here, e = {a2,N — (^t)/o'2,n is the relative difference 
between the forcing frequency of the resonant dynamic tide and the eigenfre- 
quency of the oscillation mode involved in the resonance, Q2,n is the overlap 
integral, and ^2,N is the vector of the Lagrangian displacement of the oscillation 
mode involved in the resonance. The subscript denotes the radial order of the 
mode. A generalization of Eq. Q accounting for the effects of energy exchange 
between the mass elements and their surroundings may be found in Willems, 
van Hoolst, & Smeyers (2003). 

It follows that even though the time-dependency of the displacement field 
is governed by the forcing frequency imposed by the resonant dynamic tide, the 
basic properties and general behavior of (f^) are determined by the eigen- 
function of the oscillation mode involved in the resonance. The amplitude of 
the displacement field is furthermore determined by the inverse of the relative 
frequency difference e and by the overlap integral Q2,n which is proportional to 
the work done by the tidal force through the excited mode (for a precise defini- 
tion see, for example, Zahn 1970, Press & Teukolsky 1977, Smeyers et al. 1998). 
For main-sequence stars, the behaviour of the overlap-integral is such that the 
coupling between the oscillation mode and the tidal force generally weakens with 
increasing radial order of the mode. 

The effects of resonant dynamic tides on the appearance and evolution of 
a binary may be quite substantial. Firstly, the enhanced amplitude of the tidal 
motions and the associated radial- velocity variations provide the most favourable 
circumstances for the detection of tidally induced oscillations in close binaries. 
Willems & Aerts (2002) have shown that besides the increase in amplitude, the 
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radial-velocity variations associated with a single resonance also become more 
and more sinusoidal when the forcing frequency becomes closer to the eigen- 
frequency of the mode involved in the resonance. Secondly, the enhanced tidal 
action can significantly accelerate the secular evolution of the binary's semi- 
major axis and the star's rotational angular velocity (Savonije & Papaloizou 
1983, 1984; Willems et al. 2003). Since the forcing frequencies induced in the 
star depend on n and 0,, one would therefore expect resonances to be fairly 
short-lived as the binary evolves rapidly through and away from them. Witte 
& Savonije (1999, 2001), however, have shown that when stellar and orbital 
evolution are taken into account simultaneously, the changes in the forcing fre- 
quencies and the eigenfrequencies may compensate one another and the binary 
may become locked in a resonance for a prolonged period of time. 

Despite the possibility of resonance lockings, firm evidence for the presence 
of resonantly excited oscillation modes in components of close binaries is still 
scarce. In recent years, however, systematic large-scale surveys of pulsating stars 
in close binaries containing early-type stars have yielded some very promising 
candidates in the search for tidally induced oscillations. The stars HD 177863, 
HD 209295, and HD 77581, for example, are all multi-periodic oscillators with 
observed frequencies equal to an integer multiple of the orbital frequency^. The 
binary properties and possibly resonant frequencies of these three systems are 
summarized in Tabled HD 209295 is currently believed to have either a neutron 
star or white dwarf companion (Handler et al. 2002), while the companion to 
HD 77581 is the well-known Vela X-1 pulsar (Quaintrell et al. 2003). The nature 
of the companion to HD 177863 is currently unknown. 



Table 1 . Examples of stars oscillating with observed frequencies dobs equal 
to integer multiples of the orbital frequency n. 



Name 


Ml 


M2 


^orb 


e 


CTohs/n 


Refs. 


HD 177863 


3.5M0 


I.OM0 


11.9d 


0.60 


10 


1, 2, 3 


HD 209295 


I.8M0 


1.5M0 


3.11d 


0.35 


3, 5, 7, 8, 9 


4 


HD 77581 


23-29M0 


1.8-2.4M0 


8.96 d 


0.09 


1, 4 


5 



References: (1) De Cat et al. 2000; (2) De Cat 2001, (3) Willems & Aerts 2002; (4) 
Handler et al. 2002; (5) Quaintrell et al. 2003. 



4. Apsidal motion 

In close binaries, the tidal distortion of a star by its companion gives rise to a 
perturbation of the external gravitational field which in turn causes a secular 
motion of the line of the apses. The most commonly used formula for the rate of 
secular apsidal motion was derived first by Cowling (1938) and subsequently by 



^Note that in the observer's (non-rotating) frame of reference the condition for a mode with an 
observed frequency cjobs to be resonantly excited is CTobs — fen. 
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Sterne (1939) under the assumption that the orbital period is long in comparison 
to the periods of the free oscillation modes of the component stars. The formula 
is given by 

where vj is the longitude of the periastron, k2 is the apsidal-motion constant, 
and 

/(.^) = (l-e^)-^(l4.^ + i.'). (4) 

The apsidal-motion constant k2 depends on the internal structure of the star and 
measures the extent to which mass is concentrated towards the stellar center. 
The constant takes the value ^2 = in the case of a point mass and the value 
k2 = 0.75 in the case of the equilibrium sphere with uniform mass density. For 
main-sequence stars, /c2 is typically of the order of 10~^ — 10~^. 

Within the framework of the theory of dynamic tides, the rate of secular 
apsidal motion is determined by adding the contributions to the apsidal-motion 
rate stemming from the various partial static and dynamic tides generated by the 
individual terms in the expansion of the tide-generating potential. The resulting 
rate can be cast in a form similar to Eq. @ as 

(Claret & Willems 2002). The generalized apsidal-motion constant A;2,dyn now 
renders the response of the star to the forcing frequencies induced by the tide- 
generating potential and therefore depends on the rotational angular velocity fi, 
the mean motion n, and the orbital eccentricity e. The constant is determined 
by numerical integration of the system of differential equations governing forced 
nonradial oscillations of a spherically symmetric equilibrium star. 

Unlike the classical apsidal-motion constant /c2, the generalized apsidal- 
motion constant A;2,dyn can be negative as well as positive so that periastron re- 
cessions can occur as well as periastron advances. Periastron recessions typically 
arise for close resonances when the forcing frequency of the resonant dynamic 
tide is larger than the eigenfrequency of the oscillation mode involved in the 
resonance. When the forcing frequency of the resonant dynamic tide is smaller 
than the eigenfrequency of the oscillation mode involved in the resonance, the 
apsidal motion is generally a periastron advance. 

In the limiting case of long orbital and rotational periods, the rates of 
secular apsidal motion derived by Cowling (1938) and Sterne (1939) agree with 
the formula established within the framework of the theory of dynamic tides 
up to high orbital eccentricities (Quataert, Kumar, & Ao 1996; Smeyers et al. 
1998). The agreement rests on the property that when all forcing frequencies ut 
tend to zero, the generalized apsidal-motion constant A:2,dyn tends to the classical 
apsidal-motion constant k2- For binaries with shorter orbital and/or rotational 
periods, deviations arise due to the increasing role of stellar compressibility at 
higher forcing frequencies and due to resonances of dynamic tides with free 
oscillation modes (Smeyers & Willems 2001). 
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Figure 1. Variations of the relative difference A as a function of the orbital 
period for a 5 Mq (left-hand panel) and a 20 Mq (right-hand panel) ZAMS 
star, and for the orbital eccentricities e = 0.25 and e = 0.5. The rotational 
angular velocity was set to O = 0.01 n so that the forcing frequencies are 
predominantly determined by the mean motion n. The response of the star 
to the forcing frequencies was determined in the isentropic approximation. 
The amplitudes near resonances are artificially limited to a finite value as 
explained in Smeyers & Willems (2001). 



A useful quantity to illustrate the deviations between the apsidal-motion 
rate derived by Cowling (1938) and Sterne (1939) and the rate derived within 
the framework of the theory of dynamic tides is the relative difference 

^ _ (dCT/dt)clas - {dw/dt)dyri _ k2 - fc2,dyn 

{dzu/dt)dy, 

The variations of A as a function of Porb are displayed in Fig. ^ for a 5Mq and 
a 20 Mq ZAMS star, and for 0, = 0.01 n. The relative differences are mostly 
negative so that for binaries with shorter orbital periods the classical formula 
yields somewhat too small values for the rate of secular apsidal motion, and 
thus somewhat too long apsidal-motion periods. The peaks observed at shorter 
orbital periods are caused by resonances of dynamic tides with free oscillation 
modes of the tidally distorted star. These peaks are superposed on a basic 
curve which represents the systematic deviations caused by the increasing role 
of the stellar compressibility at shorter orbital periods. 

The systematic deviations due to the stellar compressibility increase with 
increasing mass of the ZAMS star. This behaviour is related to the decrease in 
the amount of central condensation and the associated increase in the impact 
of tides with increasing stellar mass. The systematic deviations also become 
larger with increasing orbital eccentricity due to the larger number of higher- 
frequency tides contributing to the tide-generating potential. The number of 
resonant dynamic tides on the other hand tends to decrease with increasing 
stellar mass. The reason for this is that ZAMS stars with higher masses have 
smaller radiative envelopes which causes their eigenfrequencies to be larger and 
more widely spaced than those of a lower- mass ZAMS star. For a more extended 
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discussion of the relative differences A for ZAMS stars, we refer to Smeyers & 
Willems (2001). 

The extent of the deviations caused by the compressibihty of the stellar fluid 
also depends on the evolutionary stage of the star (Willems & Claret 2002) . The 
dependency is primarily through the evolution of the radius and the associated 
change of the star's dynamical time scale. In particular, as the star evolves on 
the main sequence, the radius and the dynamical time scale increase so that in 
relative terms the orbital period becomes shorter in comparison to the star's 
dynamical time scale. Correspondingly, the forcing frequencies become larger 
when expressed in units of the inverse of the star's dynamical time scale. The 
deviations due to the compressibility of the stellar fluid are therefore larger for 
a model with a larger radius. The evolution of the star furthermore also affects 
the deviations caused by the resonances of dynamic tides with free oscillation 
modes of the component stars. In particular, the effects of the resonances tend 
to be larger for stars near the end of core-hydrogen burning than for stars on the 
zero-age main sequence. This effect is again related to the size of the radiative 
envelope, which increases which increasing age of the star on the main sequence. 

The dependency of the apsidal-motion rate on the stellar structure makes 
the comparison between theoretically predicted and observationally derived ap- 
sidal-motion rates an excellent tool to probe the interior of close binary compo- 
nents (see, e.g.. Claret &; Gimenez 1993). In order to reliably compare theory 
and observations, the contributions to the apsidal-motion rate due to stellar 
rotation and general relativity must be taken into account in addition to the 
contribution due to tides. 

Stellar rotation contributes to the apsidal motion in two ways. Firstly, the 
distortion of the star due to rotation contributes to the perturbation of the 
external gravitation field. The associated rate of secular apsidal motion again 
depends on the star's internal structure through the apsidal-motion constant k2 
and on the stellar and orbital dimensions through the factor [Ri/af'. The rate 
is given by 

(e.g. Sterne 1939). Secondly, rotation tends to increase a star's central conden- 
sation and therefore decreases the apsidal-motion constant k2- Claret (1999) 
showed that the difference in log k2 can be approximated as 

The general relativistic contribution to the rate of secular apsidal motion 
on the other hand is independent of the star's radius and internal structure. 
Instead, it depends on the total system mass Mi + M2, the orbital semi-major 
axis a, and the orbital eccentricity e: 



dw \ _ 3G 2t^ M1 + M2 
dt /gr c2 Porb a(l - e2) 



(9) 



(Levi-Civita 1937, Peters 1964). The general relativistic apsidal-motion rate is 
often small, but may become dominant in binaries with high orbital eccentricities 
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Figure 2. Variations of the relative difference A as a function of the ratio 
of the rotational angular velocity to the orbital angular velocity at periastron 
for the binaries WPyx and GGOri. 



(e.g. V1143Cyg with Mi = IAMq, M2 = I.SMq, Porb = 7.6days, and e = 
0.54). 

Testing theories of stellar structure and evolution by comparing predicted 
and observed apsidal-motion rates furthermore requires both accurate obser- 
vations and high-quality stellar models (Claret & Willems 2002). Due to the 
strong dependence of the apsidal-motion rate on the ratio Ri/a, the observations 
must yield accurate determinations of the stellar and orbital dimensions, while 
the stellar models must be able to simultaneously match the observationally in- 
ferred masses, radii, and effective temperatures of both stars at the same age. 
In the context of the theory of dynamic tides, the orbital and rotational periods 
must furthermore be known with a high degree of accuracy in order to precisely 
determine the forcing frequencies induced in the binary components. While the 
accurate determination of the orbital period usually poses no problem, the de- 
termination of the rotational period of close binary components often remains a 
challenge. 

The lack of solid rotational angular velocity determinations may yield large 
uncertainties in the calculation of theoretical apsidal-motion rates within the 
framework of the theory of dynamic tides. This is illustrated in Figs. |21 and 
|51 where the relative differences A for the eclipsing binaries VVPyx, GGOri, 
UOph, and V451 Oph are displayed as a function of the rotational angular 
velocity in units of the orbital angular velocity fip at periastron. In the 
cases of VV Pyx and GG Ori, the relative differences are small and entirely due 
to the effects of the stellar compressibility for the entire range of rotational 
angular velocities considered. In the cases of UOph and V451 Oph, on the 
other hand, the answer to the question of whether or not dynamic tides play a 
role in the determination of the contribution of tides to the apsidal-motion rate 
strongly depends on the value of the star's rotational angular velocity. It can 
furthermore be seen that a small error in 0, can mean the difference between 
being in a resonance or not. 

Claret & Willems (2002) showed that when all the aforementioned uncer- 
tainties and contributions to the apsidal-motion rate are taken into account, 
a generally good agreement between theoretically derived and observationally 
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Figure 3. Variations of the relative differenee A as a function of the ratio 
of the rotational angular velocity to the orbital angular velocity at periastron 
for the binaries U Oph and V451 Oph. 



inferred apsidal-motion rates is found. The authors furthermore found no sys- 
tematic deviations between predicted and observed apsidal-motion rates as a 
function of the evolutionary stage of the stars. Despite the good agreement, one 
should keep in mind though that the theoretical treatment of the problem is far 
from complete. The detailed effects of rotation on the tidal motions and the 
internal structure of the stars, for example, has yet to be consistently included 
in the calculations. Claret &: Willems (2003) performed a first preliminary ex- 
ploration of the effects of the modified internal structure and found that huge 
differences may arise from using rotating stellar models instead of non-rotating 
ones. Their exploration, however, did not take into account the effects of the 
Coriolis force, which form an indispensable ingredient to the problem but also 
significantly complicate the theoretical treatment (see Rocca 1987, 1989). 



5. Conclusions 

In recent years, significant progress has been made in the theory of dynamic tides 
in close binaries as well as in the observations of tidally active close binaries. 
Prom the theoretical side, the effects of dynamic tides have been incorporated 
into the classical apsidal-motion test to stellar structure and evolution, while 
from the oservational side more and more evidence is arising supporting the 
theoretically predicted resonant excitation of free oscillation modes. As more 
and more data becomes available from large-scale surveys of pulsating stars in 
close binaries (e.g. Harmanec et al. 1997; Aerts et al. 1998, 2000), the full 
potential of tides as a tool to test theories of stellar and binary evolution may 
soon be within our reach. The known sample of eclipsing binaries showing 
apsidal-motion may furthermore significantly benefit from the large number of 
new systems expected to be found from large-scale surveys looking for exoplanets 
transiting the disk of their host star (see Horne 2002 and Charbonneau 2003 for 
an overview of ground-based transit searches). 

Acknowledgments. BW is grateful for the financial support supplied by 
the organizing committee and NASA ATP grant NAG5-13236 to Vicky Kalogera. 
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